Spin lattices with two-body Hamiltonians for which the ground state encodes a 

cluster state 
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We present a general procedure for constructing lattices of qubits with a Hamiltonian composed 
of nearest-neighbour two-body interactions such that the ground state encodes a cluster state. We 
give specific details for lattices in one-, two-, and three-dimensions, investigating both periodic and 
fixed boundary conditions, as well as present a proof for the applicability of this procedure to any 
graph. We determine the energy gap of these systems, which is shown to be independent of the 
size of the lattice but dependent on the type of lattice (in particular, the coordination number), 
and investigate the scaling of this gap in terms of the coupling constants of the Hamiltonian. We 
provide a comparative analysis of the different lattice types with respect to their usefulness for 
measurement-based quantum computation. 

PACS numbers: 03.67.Lx 



I. INTRODUCTION 

There is currently considerable interest in preparing 
exotic quantum states of many-body systems which can 
be used as resource states for measurement-based quan- 
tum computation (MBQC) - that is, quantum computa- 
tion that proceeds solely through local adaptive measure- 
ments on single quantum systems [11, 0, S 0, H, @| . The 
canonical example of such a resource state is the clus- 
ter state which is a universal resource for MBQC 
on suitable lattices or graphs It may be possible to 
prepare such a cluster state dynamically in atomic sys- 
tems such as an optical lattice @ or using single pho- 
tons However, one exciting possibility is that 
such resource states might be the non-degenerate ground 
state of a "natural" Hamiltonian lattice system. If the 
system is gapped, then one simply needs to cool it suf- 
ficiently in order to obtain the desired state (although, 
even f or g apped systems, this cooling process may be dif- 
ficuh ^). 

Consider the cluster state on a lattice £, defined as the 
unique -1-1 eigenstate of a set of stabilizer operators 5*^ — 
Xfj_ ^w^^ Z^, where (Z^) is the Pauli X (Z) operator 
at site fi and where v ^ fj, denotes that is connected to 
/I by a bond in the lattice C. The Hamiltonian 



fj.ec 



(1) 



with A an energy constant, has the cluster state as 
its unique ground state [^. In addition, this system 
is gapped (the gap is 2A), and such a system can be 
cooled efficiently [13] ■ However, for any non-trivial lat- 
tice or graph, this Hamiltonian involves many-body in- 
teractions, as opposed to the two-body interactions that 
occur frequently in nature. 

An obvious question, then, is whether it is possible to 
realize any given highly- entangled quantum state as the 
ground state of a Hamiltonian with only two-body in- 
teractions. Haselgrove et al. [l^] proved that this is not 



possible in general, and Nielsen [T4l| used this result to 
prove that a cluster state on a computationally universal 
(i.e., two-dimensional or higher) lattice cannot arise as 
the ground state of a Hamiltonian with only two-body 
interactions. However, investigations into quantum com- 
plexity theory [l^ [3| have demonstrated that cluster 
states (and other such states that are universal) can be 
approximated by the ground state of a local two-body 
Hamiltonian. The key idea is to make use of "mediating" 
ancilla qubits to create an effective many-body coupling 
out of two-body interactions. The problem with such 
methods is that the detailed parameters in the perturb- 
ing Hamiltonian must be controlled with a precision that 
increases with the size of the system [l7| , making such 
approaches impractical for the task of creating cluster 
states on large lattices. 

Using an alternate method based on the idea behind 
projected entangled pair states (PEPS) [18], Bartlett and 
Rudolph [l^ proved that it was possible to obtain a state 
that closely approximates an encoded cluster state on a 
square lattice using a Hamiltonian with only two-body 
nearest-neighbour interactions. In addition, they proved 
that MBQC can proceed using such an encoded resource 
state, still requiring only adaptive single-qubit measure- 
ments. 

In this paper, we present a general method for con- 
structing two-body nearest-neighbour Hamiltonian sys- 
tems for which the ground state encodes a cluster state, 
based on the techniques of [ll| . Our rigorous application 
of perturbation theory reveals errors in the calculation of 
the energy gap for the square lattice investigated in [19| 
(although these errors do not affect their key result) and 
we provide a correct treatment of this case. We also 
investigate the cluster state on a one-dimensional line 
(useful for illustration, as well as for its application as a 
quantum wire [3])j hexagonal lattice in two-dimensions 
- a universal MBQC resource with the best scaling of the 
energy gap in perturbation, and the cubic lattice in three- 
dimensions - a resource state for which fault-tolerance 
thresholds have been found [13, [2l[ . We explicitly char- 
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acterise the effects of fixed boundary conditions on the 
lattice, proving that such boundary conditions do not af- 
fect the main result. Finally, we provide an outline of a 
proof that this method yields an encoded cluster state as 
the ground state on any graph. 

II. A PEPS HAMILTONIAN 

Our general method relies on the fact that the cluster 
state is simply represented as a projected entangled pair 
state (PEPS), also known as a valence-bond solid state. 

A. The PEPS representation of a cluster state 

The PEPS representation [3| is a powerful and often 
compact method of describing the state of a many-body 
system. Consider a regular lattice £ of qubits, with co- 
ordination number c (i.e., c bonds connect each qubit to 
other sites on the lattice). A PEP state on C can be 
constructed by assigning a pair of virtual quantum sys- 
tems of dimension D to each bond on the lattice, each 
pair prepared in a maximally-entangled state, and then 
applying a projection P to the c virtual systems associ- 
ated with each site. The cluster state (and a wide variety 
of other states of interest) require only D = 2 for their 
representation; in what follows, we restrict our attention 
to this case where the virtual systems are qubits. In ad- 
dition, we choose the maximally- entangled state of these 
virtual qubits to be the two-qubit cluster state 

|C2) = ^(|0)|+) + |l)|-)), (2) 

where |±) = '^^\^) ^ I-*-))- With this convention, the 
cluster state has a simple PEPS representation [3| cor- 
responding to the projection operator 

Pl = |0l)(00...0| + |U)(11...1|, (3) 

at each site, with c zeros (ones) in (00 . . . 0| ((11 . . . 1|), 
and the states |0l) and |lz,) forming a basis for the re- 
sulting qubit at each site. 

As an example, consider the PEPS representation 
of the cluster state on a two-dimensional square lat- 
tice. There are four bonds emanating from every site 
in a square lattice, and so each site possesses four vir- 
tual qubits. Virtual qubits connected by a bond are 
placed in the state IC2), and then a projection Pl = 
«)sitcs(|OL) (00001 -I- |1l)(1111|) is applied. The resultant 
state, \4>) = ^'L(8ibonds|C2) is a cluster state on the square 
lattice. 



B. A two-body PEPS Hamiltonian 

The essential idea of the method presented in this pa- 
per is to mimic the PEPS construction procedure using 



a physical two-body Hamiltonian, wherein the "virtual" 
qubits are real physical systems and the resulting PEPS 
state is encoded into logical qubits. Consider a regular 
lattice. Let C denote the set of sites, each with coordina- 
tion number c. We assign c qubits to each site, and label 
with a double index (/i, «), with fi Cz C and i = 1, 2, . . . , c. 
(The choice of the second label i is completely arbitrary.) 
Let cr^^ and a^^ denote the Pauli X and Z operators 
for the ith qubit at site /i G £. Following the PEPS 
construction, if a site fj, is connected to a site by a 
bond (denoted ^ ^ u), then we associate qubit (/j,, i) and 
{ly, i) for some i to this bond. (We note that this notation 
can become problematic with certain periodic conditions, 
but it should be clear from the context how to adjust it 
appropriately.) 

Our PEPS Hamiltonian is defined as follows. At each 
site, we require a two-body Hamiltonian with a two- 
dimensional ground state space spanned by |00 . . . 0) and 
|ll...l). For this, we choose a site Hamiltonian Hq 
which is of Ising form 

tj.ec i^j 

where i j denotes that qubits i and j are connected 
according to some graph structure. Aside from being con- 
nected, the specific form of this graph is relatively unim- 
portant; however its structure will affect the energy lev- 
els of the Hamiltonian. For example, for two-dimensional 
lattices, it is natural to choose a ring structure. 

Between sites, we define a different two-body interac- 
tion of the form 

c 

^^-T.T.'^M^'^lm)' (5) 

fiGC i=l 

where i/(z) is the site connected to fi via bond i. With 
this Hamiltonian, on every bond in the lattice ^ v 
there are two terms: cr'f ® erf x and erf ® a'f x. 
Note that the terms in V stabilise IC2), and therefore 
®bonds|C'2) IS the grouud state of V. This product of 
maximally- entangled states is the starting point of the 
PEPS construction. The site Hamiltonian Hq is meant 
to "implement" the PEPS projection by ensuring that 
the qubits at a site act as a single logical qubit; to do so, 
the site Hamiltonian Hq must be much stronger than the 
bond Hamiltonian V . One is then lead to consider the 
ground state of the Hamiltonian 

H = gHo + XV , (6) 

where g ^ X > 0, which is suitable for perturbative 
analysis in X/g. 

In jl9j, this procedure was applied to a square lattice. 
The terms in the perturbation combine at higher orders 
to yield the stabilisers of the logical cluster state, and 
the resulting low-energy theory of the lattice is governed 
by an effective Hamiltonian of the form of Eq. ^ . Fur- 
thermore, the gap to the next excited state is finite and 
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independent of the size of the lattice. Of course, because 
this is a perturbative approach there will now be correc- 
tions to the unperturbed logical eigenstates which will 
not be in the logical space. So the exact cluster state 
will not be obtained. However, these errors will be small 
(occurring with probability (A/.g)^, as we will show) and 
so a state arbitrarily close to the cluster state can be ob- 
tained. In what follows, this procedure is generalised to 
other lattice structures important in quantum computa- 
tion. 



C. General properties of the PEPS Hamiltonian 

1. Duality transformation to uncoupled sites 

We now present a simple duality transformation that 
maps the Hamiltonian ([6]) to one describing uncoupled 
sites. Consider decomposing (|6]) as a sum of commuting 
terms 77^, as 



where 



(7) 



(8) 



Note that [H^^Hi,] — for ^ ^ v. Define the unitary 
transformation CS£ to be the application of a CSIGN 
gate 



CSIGN 




a ^ (7 

(g) 

I® a'' 



(9) 



to every bond on the lattice. The action of this transfor- 
mation on the bond terms in the above Hamiltonian is 



(CS^)^^^^,)® 
term thus yields 



(CSr) 



Transforming each 



H' = {GSc)H,{CSc) 



-sE^ 



(10) 



i=l 



which is localized entirely to the site /i. Thus, the duality 
transformation CS£ yields a Hamiltonian of uncoupled 
sites, where each site Hamiltonian takes the form of a 
transverse-field Ising model on some connected graph. 

With this mapping, the spectrum of the Hamilto- 
nian ^ can be calculated explicitly, with the Hamilto- 
nian term H'^ at each site, for example by using a Jordan- 
Wigner transformation. We note at this point that each 
Hamiltonian H'^ has a nondegenerate ground state for 
all A > 0; thus, our PEPS Hamiltonian on the full lattice 
will also possess a nondegenerate ground state for A > 0. 



2. Encoded stabilizers: Constants of motion 



For each site ^, define the operator 



(11) 



1=1 



That is, is the tensor product of cr^ for every qubit at 
site together with on every neighbouring site. It is 
straightforward to show that all such operators commute 
with the Hamiltonian of Eq. ([5]), 



and with each other. 



(12) 



(13) 



Thus, if H has a nondegenerate ground state (as is the 
case for the PEPS Hamiltonian with A > 0), it must also 
be a simultaneous eigenstate of all operators K^. 

Using the duality transformation CS£ defined above, 
we find that 



(CSc)i^4CSc) = (>?)(T 



ip-.i) 



(14) 



Using the well-known solution to the transverse-field 
Ising model with Hamiltonian (fTU|) . we find that the 
ground state for A > is the +1 eigenstate of this oper- 
ator. (This can be inferred by the fact that the ground 
state is clearly the -|-1 eigenstate of (8)i=icr^^ in the hmit 
\/g oo.) Thus, we have that the ground state of the 
PEPS Hamiltonian for A > is the simultaneous -1-1 
eigenstate of all operators K^, ^ e £. 

The operators K^, then, can be viewed as encoded clus- 
ter stabilizers, and the ground state for A > as an en- 
coded cluster state. Unfortunately, for A > 0, this encod- 
ing is no longer in the ground state space of Hq spanned 
locally at sites by the states |00 . . . 0) and |11 . . . 1). Our 
perspective is to consider the encoding to be fixed in this 
space and view the ground state of the PEPS Hamil- 
tonian as the desired locally-encoded cluster state plus 
perturbative corrections in X/g. These concepts are best 
illustrated through a simple example. 



III. EXAMPLE: A 1-D LINE 

First we illustrate this approach on the simplest lattice: 
a one-dimensional line of qubits with periodic boundary 
conditions. We demonstrate explicitly that the pertur- 
bative procedure yields a low-energy effective Hamilto- 
nian on the logical qubits of the form of Eq. and 
that an approximate cluster state is obtained as the non- 
degenerate ground state. The basic steps outlined in this 
example for the one-dimensional line, suitably general- 
ized, will be applicable to more complex lattice struc- 
tures. 



■■■■a--o-a--o--Q--o- 



Insert 2 physical qubits 
at every lattice site: 



lattice site (a logical qubit) 



physical qubits 



FIG. 1: The logical lattice and the physical lattice structure 
for the ID line. 



Consider a one-dimensional line consisting of Ng qubit 
sites with periodic boundary conditions (i.e., a ring). The 
coordination number of this lattice is c = 2, and thus 
our construction requires two qubits to be placed at each 
site. This lattice structure is illustrated in Fig. [TJ The 
Hamiltonian is that of Eq. 



A. The unperturbed spectrum 

We first investigate the energy eigenvalue spectrum of 
the unperturbed Hamiltonian gHo. Because gHo is a sum 
of terms, each of the form (g) cr^ acting on a single site, 
the energy spectrum can be determined by analysing each 
site individually. At a single site, there are two energy 
levels. The ground-state is degenerate, two-dimensional, 
and spanned by the states 



100) 



10, 



111) 



(15) 



The ground state space of the unperturbed Hamiltonian 
at each site, then, can be viewed as a logical qubit. Note 
also that this ground state space is, by construction, the 
logical subspace for the cluster state PEPS projection. 
The energy of this ground state space is —g. The first 
excited state at each site is also two-dimensional, has an 
energy of g, and is spanned by the states |01) and |10). 

With the spectrum of gHg at each site, we now describe 
the spectrum of the unperturbed Hamiltonian gHo on the 
entire lattice. The lattice ground-state space is spanned 
by product states of all of the individual sites in the 
ground state (i.e. in the logical space). This ground-state 
space has energy Eq^^ = —gNs, is 2^^-dimensional, and 
is spanned by all logical states of Ng qubits. We denote 
this space Hl- The first-excited space is {2Ns • 2^^~^)- 

dimensional, and has energy e[^^ = —g{Ns — 2). Thus, 
for the unperturbed Hamiltonian gHo, the gap from the 
ground to first-excited space is 2g. The second-excited 

space has energy £^2°^ = —gi^s — 4), and so on. These 
energies will serve as the zeroth-order energies in pertur- 
bation theory for the total Hamiltonian. 



B. Perturbation theory 




Biergy 



+g 



FIG. 2: The effect of the terms in V on a single site. 



show that this term lifts degeneracy of the ground state, 
and that the logical cluster state arises as the unique 
ground state (although we also show that there are per- 
turbative corrections to this state). For details of our use 
of perturbation theory and notation, see the Appendix. 
Let the nth-order energy correction to the jth state in 

(n) 

Hl^s denoted by X^Eqj' . Let Pl be the projection onto 
the degenerate ground state space of the unperturbed 
Hamiltonian gHo, i.e., onto the "logical" space Hl- De- 
fine Pl ■= I — Pl to be the projection onto the "illogical 
space" (denoted Hi) and let the projection onto the jth 
unperturbed excited level be denoted Pj . 

To obtain a conceptual view of the perturbation it is 
useful to see the effect of F on a single site. The Hamilto- 
nian is a sum of terms of the form —a^ (g) a^; however, 
each of the Pauli operators in such a term act on differ- 
ent sites, and so we must consider the action of and 
cr^ separately. Because the logical space Hl spanned 
by 1 00) and |11), the action of will move a site out of 
the logical space; the action of will not, and simply 
induce a phase. The possible actions by the cr^ part of 
y at a single site are shown in Fig. [2l 

The first-order corrections to the energy are governed 
by the operator (see Eq. (|A15|) in the Appendix) 



PlVPl 



(16) 



Specifically, the first order energy corrections E^yj^ to the 
ground state are the eigenvalues of this operator. Because 
all of the terms in V contain a single , they all move 
a state in the logical space to the first excited state (i.e., 
VPl = PiVPl). Thus, PlVPl = 0, and there is no 
first-order correction to the energies. 

The second-order corrections are governed by the op- 
erator 



(2) = Pim(4°) 

PlVPiVPl 



gHor^pLVPL 



(4°) -i?r) 



(0)^ 



(17) 



We now turn to perturbation theory and determine the 
effect of the term W in the Hamiltonian We will 



where the expression has been simplified using VPl = 
PiVPl- The operator PlVPiVPl maps states from the 



5 



ground state space to the first excited space and then 
back to ground state space. By investigating the differ- 
ent ways of returning to the logical space after just two 
CT^ operations, it is clear that there are two possible con- 
tributions to this term: 

1. If cr^ (8) (T^ in y is applied twice to the same bond, 
it yields the identity. The first a'^ can be applied to 
any of the qubits and then must be applied again 
to the same qubit, so there are 2Ns of these terms. 

2. If (T^ is applied to each of the two qubits at a site 
(and corresponding a" operations to qubits in the 
neighboring sites), then the lattice remains in the 
ground state. We can apply the first cr^ to either of 
the two qubits at the site and so there will be two 
of these terms that occur at each site. Explicitly, 
this case will be a term applied to the logical space 
of the form 



(18) 



1=1,2 



where fi v. 

The operator 5^ := PlK^Pl^ which acts only on 
the logical space, can be determined explicitly as 
follows. Note that the product of two cr^ opera- 
tors on a single site /i (one on each physical qubit), 
restricted to the logical space, is equivalent to a 
logical X operator 



(19) 



Also, a single operator acting on either of the 
two physical qubits at a site restricted to the 
logical space, is equivalent to a logical Z operator 



(20) 



Thus, S^j, = X^j^®i,^^j^Zy. This operator is a (logical) 
stabilizer of the cluster state on this lattice. 



Therefore, we have 



PlVPiVPl = 2NsPl 



2^5,, 

M6£ 



(21) 



Substituting this result into Eq. p?| and using 



(0) 



E' 



(0) 



~2g gives 



'iNsPL 



-25 



The energies E^^^ are the eigenvalues of 6*^^^ and the cor- 
responding eigenstates of 0*^^^ will be the zeroth-order 
energy eigenstates after the degeneracy is lifted. 

Next, we identify the basis which diagonalises 0*^^-'; 
this is straightforward given the expression (P^ . As 
the cluster state is the simultaneous -1-1 eigenstate of all 



stabilizer operators 5'^, the logical cluster state on this 
lattice, denoted by |C), is an eigenstate of 0^^^. Simi- 
larly, the other eigenstates of 9^"^^ are also just the si- 
multaneous eigenstates of the stabilizers (as all such 
stabilizers commute). Explicitly, let \C{a, f3, . . .}) de- 
note the logical cluster state with a logical Z-operator 
(called a Z -error) applied to the sites a, (3, ... & C. Us- 
ing the anti-commutation relations of the Pauli matri- 
ces, \C{a,l3, . . .}) is the —1 eigenstate of SajSp, ... and 
the -1-1 eigenstate of S*^ for fi ^ a, (3,.... Therefore, 
the 2-^^ states of the form \C{a,(3, . . .}) will be eigen- 
states of 0^^-' . Furthermore, these states are orthogonal, 
as each pair of states will have a differing eigenvalue for 
least one of the 5^ operators. In summary, the set of 
states {|C), |C{a}), \C{a,f3}), . . .}, running over logical 
Z-errors at all possible sites, forms an orthogonal basis 
of Hl and diagonalises 9'^^\ 

The eigenvalue spectrum of 6^^^ is then straightforward 
to calculate using the properties of stabilisers. From the 
form of O^^") in Eq. (P^ . the lowest energy eigenstate will 
be the cluster state |C), because it is an eigenstate of 
all stabilisers in the sum J2fj.£C "^m with eigenvalue -1-1. 
Thus the second-order correction for the energy associ- 
ated with the cluster state is 



2p(2) _ 2Ns 



2Ns 



X2 

A2 ^ ^2Ns — 
9 



(23) 



Next, consider a state |C{a}) Zq,|C), a cluster state 
with a single Z-error at the site a. This state is also 
an eigenstate of all stabilizers in the sum with 
eigenvalue -f 1 except the stabilizer Sa for which it has 
eigenvalue —1. Therefore, 



X2p(2) 

^ ^|C{a}> 



2Ns + 2{Ns - 2) 
-2.9 



A2 ^ -2{Ns - 1)- 
9 



(24) 

Because there are Ns states of the form |C{q!}), this 
1*^* excited space is A^g-dimensional. Similarly, the nth 
excited space up to n — Ns is (^^) -dimensional and (to 
zeroth order) is spanned by states obtained from |C) by 
n logical Z-errors. 

Higher order corrections can be calculated by following 
a similar procedure. As noted in Sec. Ill C 11 this Hamil- 
tonian can be easily solved exactly, with a ground state 
energy given by 



(25) 



(22) There IS an energy gap 




2Xyg + 0{X'/g'), 



(26) 



to the first excited space; all higher levels have energy 
En — E^c) +nA. Note that A is independent of Ns, the 
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size of the lattice. Intuitively, then, one may associate 
logical Z errors on any site with a fixed energy A each. 

In summary, we have shown that the non-degenerate 
ground state of the Hamiltonian H = gH^ + is the 
cluster state, to zeroth order in A/(/, with an energy gap 
to the second excited state scaling as ^ A^/g. 



C. Perturbative corrections to the ground state 

We have shown that, to zeroth order in A/g, the ground 
state of the system is the logical cluster state |C). How- 
ever, the perturbation will also modify the energy eigen- 
states from their unperturbed states. To first order in V ^ 
the perturbed ground state |£'o) is given (up to normal- 
ization) as 



|i?o)(x|C)+A( 



m\c) 



E 



where 



\l)(iHi^.\l)^\C) ^|C) 



(/|g(3)|c) 



(2) 



lO), (27) 



(28) 



For this perturbation, fl^^) = Q; however, there exist 
states \j) £ Hi such that {j\V\C) ^ 0. 

Note that is a sum of 2Ns terms of the form (g) 
acting across a bond. Each of these terms applied to |C) 
gives an excited state of the form 



(29) 



Using the anti-commutation relations of the Pauli matri- 
ces, the terms in Hq act on |fc(^t,i)) as 



(30) 
(31) 



Hence gHo\k^^,,)) = -g(iV5 - 2)|fc(^,,)) = ^ 



and therefore the states \k 



(p..i)l 



are in the first excited 



space of gH^. Eq. (PT|) also shows that and 
1^(^2,22)) Ml 7^ are eigenvectors with different 
eigenvalues for the operator a^^^ ® a^^^ 2) ^^'^ thus 
they are orthogonal. However, recalling from earlier 

that Kf^ := (Eif^icr^^ (8) crf^(j) i) stabilises |C), we have 
that (fc(^,i)|fc(^,2)) = {C\Kf,\C) = 1 and thus |fc(p,i)) = 
|A(/x,2))- Hence |F|C) = 2, and 



l^^o) cx |C) + A 



E 



p(o) 

E l^(M,l)> 



(0) '^(A-a) 



(32) 



There are iVg states in the above sum, which determines 
the normalization. Thus, we can calculate the fidelity 
F = |(C|i?o)P of the ground state with the exact cluster 
state, which in this case is found to be 



F 



1 



l + NsXyg^ 



(33) 



This fidelity decays rapidly for increasing Ns, which is 
unsurprising given that it is comparing quantum states 
on a large lattice and is an extensive quantity. For any 
lattice system with Ns large, this fidelity is known to 
scale as = , where d is an intensive quantity that 
can be interpreted as the average fidelity per site [23j . 
Precisely, 



logd :— lim 

Ns^o 



NzHogF, 



(34) 



which is found to satisfy 

logd=-iV5ilog(l + 7VsAV5') 
>-iV^ilog(l + AV5T" 
= -log(l + AV5')- 



(35) 



Thus d > (1 -I- A'^/g^)"^, which is independent of Ns- 
This result demonstrates that the ground state is "close" 
to the ideal cluster state, as quantified by a large average 
fidelity per site, for X g. 



IV. UNIVERSAL RESOURCES FOR MBQC 

Although it serves as an illustrative example of the 
techniques presented in this paper, the cluster state on 
a line is not a universal resource for MBQC; a higher- 
dimensional lattice is required. In this section, we apply 
the perturbative procedure to lattice structures that are 
interesting from a MBQC perspective, and comment on 
their utility. 



A. Hexagonal lattice 

For a hexagonal lattice in two dimensions with Ns sites 
and periodic boundary conditions, the coordination num- 
ber is 3, and we require three physical qubits at each site 
(see Fig. The Hamiltonian for the lattice is again 
given by Eq. 

We investigate the spectrum of the unperturbed 
Hamiltonian gHg by considering its action at a single 
site, where the three site qubits interact via the Ising 
coupling on a ring. The ground-state is degenerate, two- 
dimensional, and spanned by the states 



|000) |0, 



1111) 



(36) 



which encode our logical qubit. The energy of this 
ground state space is —3g. The first excited state is 





Biergy 

+g 





-3g 



FIG. 3: The hexagonal lattice structure, and the effect of the terms in V on a single site. 



six-dimensional, and has an energy of g. Thus, for 
the entire lattice of Ns sites, the ground-state space 
has energy Eq^^ = —3gNs, is 2^^ -dimensional, and is 
spanned by all logical states of Ns qubits. The first- 
excited space is {6Ns ■ 2^^~^)-dimensional, and has en- 
ergy Ef"^ ^ ~g{3Ns - 4.). 

We now turn to perturbation theory. It is again useful 
to obtain a conceptual view of the effect of y on a single 
site, as illustrated in Fig. [31 As the ground state space is 
spanned by |000) and at each site, only the action of 
tr^ (and not cr^) will map states out of the logical ground 
state space. The possible actions by the part of V at 
a single site are shown in Fig. [31 Once again, PlVPl = 
and there is no first-order correction to the energies. 

The second order corrections E, ' 
of the operator 
operator PlVPiVPl, we examine Fig. [31 and the ways of 
returning to the logical space after just two applications 
of V. It is clear that there is only one possible contribu- 
tion: if a tJ^ (8) cr^ in V is applied twice to the same bond, 
this will yield the identity. The first tr^ can be applied 
to any of the qubits and then must be applied again to 
the same qubit, so there are 3Ns such terms. Hence 



are the eigenvalues 
2' defined in Eq. (jXTT)) . To evaluate the 



PlVPiVPl-^SNsPl 



(37) 



Using this result in Eq. ^ as well as E^l^^ - e'^^ = -Ag 
gives 



to the ground-state energy - an energy shift - given by 



\'eP 



4g 



(39) 



The energy degeneracy of the ground state has still not 
been broken at second order and we must proceed to 
third order. 

(3) 

The third order corrections Eqj are the eigenvalues of 
the operator 9^^^ given by 



0(3)^P^y[(£;(")_i/„)-ip^y]Vi 
PlVPiVPiVPl 



(0) 



(40) 



where the expression has been simplified using VPl — 
PiVPl- With three applications of the perturbation 
V, the operator PlVPiVPiVPl maps states out of the 
ground state space and then back again via the first ex- 
cited space. Again investigating Fig. [31 it is only possible 
for the lattice to remain in the ground state after three 
perturbation terms if operators are applied to each 
of the three qubits at a site (and, through V, the corre- 
sponding operators to the qubits on each of the neigh- 
boring sites). That is, this case will be a term applied to 
the logical space of the form 



(41) 



=1,2,3 



(2) _ 3NsPl 
(-45) 



(38) 



This operator simply acts as the identity on the logical 
space and so there is a constant second-order correction 



where j/(i) is the site connected to (/i, i) by a bond. Just 
as in the case of the line, the operator acts on the 
logical space as S"^ PlK^Pl = Yliyr^^ "^^^ ^ logical 
cluster-state stabilizer operator. The three qubits at the 
site can be ordered in 3! possible ways, and so there will 
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be 3! of these terms that occur at each site. Therefore, 

(42) 



PlVPiVPiVPl = -3! ^ , 



and 



Once again, as in the hne, the set of 2^^ states 
{\C),\C{a}),\C{a, (3}), . . .}, running over logical Z- 
errors on the cluster state |C) at all possible sites, forms 
an orthogonal basis of TCl and diagonalises 9^^\ The 
cluster state |C) is the unique lowest eigenstate of 9^^\ 
The third-order correction for the energy associated with 
this state is 



r 



(44) 



Again, this case is simple enough to analyze analytically; 
the ground state of the Hamiltonian H = gHo + XV has 
energy 



\c) 



/A /A2 A 
gNs{l + - + 2J— + - 



A 
9 



1 A2 1 A3 
-3gNs(l + -— + - — 
4g^ 8g^^ 



(45) 



The nth excited space up to n = Ns is (^•^) -dimensional 
and is spanned (to zeroth order) by states obtained from 
|C) with n logical Z errors. These states have energy 
En = E]^(j^ + 71 A where 



A 



A 



^5 
3 A3 
45 




2+0(AY9') 



(46) 



We can also calculate the first-order corrections to 
the ground state |C), by finding states \j) e Tii such 
that {j\V\C) ^ 0. As before, define := cr|?^ ® 

i) 1^) ■ determining the effect of each of the terms 
in gHo on |fc(p,i)) it is clear that they are in the first ex- 
cited space of gHo and are orthogonal to each other. To 
first order in X/g, 



(47) 



Comparing this ground state with the ideal cluster state, 
we find that the average fidelity per site d is bounded by 
d>(l + 3AV(43)2)-i. 



B. Square lattice 

We now repeat the above procedure for a 2D square 
lattice with Ns sites and periodic boundary conditions. 
This case was originally examined in [l^; however, our 
detailed derivation reveals some errors in their calcula- 
tion of the perturbed energies and the gap. 

The coordination number of this lattice is 4, and so 
four physical qubits are necessary at each site (see Fig.d]). 
The Hamiltonian for the lattice is again given by Eq. ([6]), 
again with a ring of four qubits coupled via an Ising in- 
teraction. There are now three energy levels of gHo at 
a single site. The ground-state space of gHo at a single 
site is spanned by the states 



10000) =: \0l) 



mil) 11, 



(48) 



and the energy of this ground state space is —4(7. The 
first excited state is twelve-dimensional, and has an en- 
ergy of 0. The second excited state is two-dimensional 
and has a energy of 4g. So, for the entire lattice of Ns 



sites, the ground-state space has energy E^ 



(0) 



-4gNs, 



is 2^®-dimensional, and is spanned by all logical states 
of Ns qubits. The first-excited space has energy e[^'' — 
~4:g{Ns — 1) and the second-excited space has energy 

= -4g{Ns - 2). 
The possible actions by the part of y at a single 
site are shown in Fig. |4l We now follow the identical 
procedure as done previously, and find 



0(1) =0, 

^ msPL 
(-4g) 
0(3) =0, 

0(4) 5 



E^. 



16g3 ^ 



Ns 
16g' 



■Pl 



(49) 
(50) 
(51) 
(52) 



That is, there are no first- or third-order corrections 
to the energy; at second-order there is a constant en- 
ergy shift A^i^Q^) = —NsX^/g to the ground state; at 
fourth-order the degeneracy is broken. In the expres- 
sion for 6^^\ the first term is recognized as propor- 
tional to the cluster Hamiltonian: the sum of stabilis- 
ers of the cluster state S^. Therefore the set of 2^^ 
states {|C), |C{q;}), |C{q;,/3}), . . .}, running over logical 
Z-errors at all possible sites on the cluster state |C), is 
an orthogonal basis of TYl which diagonalises 9^'^^ . 

The cluster state |C) is the unique lowest eigenstate 
of 9^^\ because it is an eigenstate of all stabilisers in 
the sum with eigenvalue -1-1. The fourth-order 

correction for the energy associated with this state is 



^ ^\C) 



3.. 

8 'g^ 



(53) 



We note that this result differs, by numerical factors, 
from the result of [T9l|. (The error in [l^ arises from 




missing contributions to the perturbation operator 0^^-' 
in Eq. ([5^ .) Higher-order corrections follow in a similar 
fashion, and a complete analytic solution for the ground 
state energy is found to be 




(54) 



The nth excited space up to n = Ns is ( 
dimensional and is spanned (to zeroth order) by states 
obtained from |C) by n logical Z errors. These states 
have energy En — E\(j^ + nA, where 



A ;= 




Once again we can also calculate the first-order correc- 
tions to the ground state |C) is calculated to be 



(56) 



where Ifc, 



Comparing this 

ground state with the ideal cluster state, we again find 
that the average fidelity per site d is bounded by d > 
(l + 4AV(4.g)^)-^ 



C. Cubic lattice 

We apply the now familiar procedure to a cubic lattice 
in three dimensions with Ns sites and periodic boundary 



conditions. The coordination number is 6, and so six 
physical qubits are necessary at each site (see Fig. [5]). 
The Hamiltonian for the lattice is given in Eq. ([6]), where 
we arrange the 6 qubits on the vertices of a octahedron, 
and place (8) couplings between all qubits connected 
by an edge of the octahedron, as in Fig. [S] There are four 
energy levels of gHo at a single site; the ground-state is 
degenerate, two-dimensional, and spanned by the states 



1000000) |0i) , 1111111) =:|1l) 



(57) 



The energy of this ground state space is —12g; the first 
excited state space has energy —ig; the second excited 
state space has energy 0; the third excited state space 
has energy Ag. 

The possible actions of cr^ at a single site are shown in 
Fig. [5l Again following our general perturbative proce- 
dure we find 61(1) = 9'-^^ = 9'-^^ = and 



0(6) = 



6Ns 



Ns 
256g3 
_13Ns_ 
'491525' 



Pl, 



-Pl 



83 



16384g' 



(58) 
(59) 
(60) 



Thus, there are two constant energy shifts at second and 
fourth order of 



A2i?(2) = - 



256g3 



(61) 



In the expression for 0(6) , the de generacy is broken by the 
^ which are the cluster stabiliz- 



terms 5^ = A^ 



ers. The set of 2^^ states {\C),\C{a}),\C{a, f3}), . . .}, 
running over logical Z-errors at all possible sites of the 
cluster state |C), forms an orthogonal basis of Hl which 
diagonalises 0(6) . 
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{|100001>,...} 



{|110001>,...} 



{|011110>,...} 



V X 




{|110000>, 


X 


{|111000>,...>} 


X 


{|001111>, 






|101000>,...} 




<— » 


10101 11 >,...} 



Biergy 
4g 



V / 



{|100000>,...} 



{ 101111 1>,...} 





-4g 



-12g 



FIG. 5: The cubic lattice structure, and the effect of the terms in 1/ on a single site. 



The cluster state |C) is the unique lowest eigenstate of 
O^^^f . The sixth-order correction for the energy associated 
with this state is 



131 A'^ 

Nq — 

24576 g5 



(62) 



Therefore, to sixth order, the energy of the ground state 
is 



1 A2 



1 A4 



131 A^ 



3 • 210 5^ ' 32 . 215 g6 



(63) 

The nth excited space up to n = Ns is (^•^) -dimensional 
and is spanned (to zeroth order) by states obtained from 
|C) by n logical Z errors. These states have energy En = 



E, 



\c) 



nA, where 



A := 



83 A^ 
8192^ 



(64) 



Once again we can also calculate the first-order correc- 
tions to the ground state |C) is calculated to be 



^ nec i=i 



(65) 



where Ifc, 



Comparing this 

ground state with the ideal cluster state, we find that 
the average fidelity per site d is bounded by c? > (1 + 
6AV(8.g)2)-i. 



D. Implications for MBQC 



are all universal resources for quantum computation. In 
each case, it has been shown above that the perturba- 
tive procedure produces a non-degenerate ground state 
which approximates an encoded cluster state on the lat- 
tice. We chose to investigate each of these lattice struc- 
tures because each has a unique relevance to the study of 
MBQC. The 2D square lattice is the canonical example 
for use in cluster-state quantum computing and was the 
original lattice structure presented in [ij. This lattice is 
also the most easily accessible to experimental investi- 
gation in cold atomic systems A hexagonal lattice 
was also examined above because (as argued in [igj) the 
perturbative procedure produces a cluster state with the 
largest energy gap for a given ratio X/g. We discuss the 
implications of this observation below. Finally, recent 
work [2O, [2l| has shown that fault-tolerant thresholds can 
be found for MBQC if the lattice used is 3-dimensional. 

Following on from the discussion in [T9| . we now com- 
pare the results for each lattice and relate it to its useful- 
ness for quantum computation. There are two sources of 
error when using the ground state obtained in the pertur- 
bative procedure for cluster-state quantum computation. 

First, note that errors will arise because the ground 
state of the system is not exactly the cluster state, 
but contains perturbative corrections (cf. Eq. ([47]) . ((56)) . 

In each case the ground state is given by a su- 
perposition of the cluster state with other first-excited 
states corresponding to "errors" (E) applied to all 
bonds on the cluster state independently. This error rate 
is quantified by the average fidelity per lattice site c?, de- 
fined by Eq. which was explicitly bounded in all of 
the above examples. This bound takes the general form 



The cluster states on the three lattice types examined 
in this section (the square, hexagonal and cubic lattices) 



(66) 
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where fc is a constant of order one which depends on the 
lattice. This bound tells us that, for A ^ 5, the ground 
state is very close to the cluster state, and that the er- 
ror rate for the independent cr^ (g) is less than kX^ / g^. 
Because we require X <^ g, this error probability will 
be small. The effect, and possible error correction, for 
such two-qubit correlated errors has not yet been inves- 
tigated, but the independence and locality of the errors 
makes them amenable to existing error correction tech- 
niques. In particular, we note that such an error can be 
identified by checking sites, each of which should be in 
the code space spanned by |00 . . . 0) and |11 . . . 1). Errors 
of the form (g) tr^ will cause a correctable error to this 
code space (which must also include a phase correction 
to the appropriate neighbouring site) provided that the 
lattice has coordination number c > 1, i.e., for lattices 
of higher dimension than the 1-D line. This correction 
scheme would require measurements of multiple qubits, 
and it would be worthwhile to investigate whether such 
error correction could be performed using single-qubit 
measurements. 

The main difference arising in the calculations for each 
lattice structure, however, is the order in the perturba- 
tion theory at which the ground-state degeneracy is bro- 
ken. This occurs at third order for the hexagonal lattice, 
fourth order for the square lattice and sixth order for the 
cubic lattice. In general, the order at which perturba- 
tion theory breaks the ground-state degeneracy is given 
by the coordination number of the lattice. This result 
leads directly to a dependence of the energy gap A on 
the coordination number c of the lattice, as 



A^iX/gy. (67) 



In all cases the energy gap A is independent of the size 
of the lattice, i.e. the system is gapped. Given that the 
rate at which the thermal state of this system will ex- 
hibit Z-errors depends explicitly on the size of this gap, 
the system will be less sensitive to these errors if the en- 
ergy gap A is made larger. The hexagonal lattice will 
have the largest energy gap, as is consequently less sensi- 
tive to thermal errors. It should be noted, however, that 
methods to identify and correct for such thermal errors 
(and Pauli errors in general) within the MBQC paradigm 
currently exist only for 3-dimensional lattices [20I [2l| . 
(See also l25|].) The 2-D lattices (hexagonal and square) 
may not allow for error correction of such thermal errors 
using only single-qubit measurements; this remains a key 
open question. 

We note that there exists a trade-off between these 
two types of errors when using the state for MBQC. In- 
creasing the value of X/g will reduce the probability of 
thermal errors at a given temperature but also perturb 
the ground state away from the cluster state. 



V. FIXED BOUNDARY CONDITIONS 

The perturbative approach has so far been successful 
in producing the cluster state on each lattice type with 
periodic boundary conditions. We now analyse the effect 
of placing fixed boundaries on the lattice. 

A. A line with fixed boundaries 

We first examine a line with fixed endpoints. The in- 
terior sites still have coordination number 2, and so we 
require two physical qubits at these sites. However, the 
boundary sites will consist of just a single physical qubit. 
Denote the number of interior sites by Ns, so that there 
are (Ng + 2) sites in the entire line. In addition, denote 
the two boundary sites by the labels fi = Bi and fi = B2. 

The Hamiltonian for this lattice will remain that of 
Eq. ([5]), where we do not place any site Hamiltonian 
term on the boundary sites. The unperturbed energy 
spectrum at each of the interior sites is unchanged from 
the periodic boundary case (as in Fig. [2]). The two 
boundary qubits, however, have zero unperturbed en- 
ergy. The spectrum of the unperturbed Hamiltonian for 
the entire line is therefore a 2^^+^-dimensional ground- 
state space with energy Eq^"^ = —gNs, and is spanned 
by all logical states of {Ns + 2) qubits. The first- 
excited space is {2Ns ■ 2''^^+-'^)-dimensional, and has en- 
ergy e[°^ =-giNs-2). 

We now turn to perturbation theory. Note that, for 
the two boundary qubits, a single application of cr^ maps 
the logical space onto itself. Thus, due to the contribu- 
tions from the boundary qubits, there is now a first order 
correction to the energy 

0(1) = p^Yp^ 

= -Pl{Kb,+Kb,)Pl 

= -{Sb,+Sb,), (68) 

where 

Kb^ = fT^, ® ^1=1,1) ' (69) 

KB,=<yl,®(Tl,=Ns.2)^ (70) 

and as usual S'p = (gj/^^ Z^. In particular, the cluster 
stabilizers for the end sites are given by the product of an 
X operator on the boundary site with a single Z operator 
on its sole neighbour. 

The first-order corrections to the ground-state en- 
ergy, i?Q^\ are the eigenvalues of 6''^^^. This operator 
is diagonal in our familiar basis for Hl oi the 2^^+^ 
states {|C), |C{a}), \C{a, (3}), . . .}, running over logical 
Z-errors at all possible sites. The states in this basis 
which are the -|-l-eigenstates of both Sbi and Sb2 will 
be the lowest eigenvalues of 0^^-'. By the property of sta- 
bilisers, the subspace T g Hl which is stabilised by Sbi 
and Sb2 is 2^^-fold degenerate, and so the degeneracy 
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is reduced by a factor of 4 at first order. Tiie cluster 
state |C) is contained in this subspace, as are all states 
with logical Z-errors anywhere except on the boundary. 
Thus the lowest energy space T has a first-order correc- 
tion given by: 



XEj^ ^ — — 2 A . 



(71) 



The next highest energy level includes states which have 
a Z-error at either of the boundary sites Bi or B2 but not 
both. These 2^s+i states have AE'^^^ = 0, which is a gap 
of 2A above the space T. The second highest energy level 
will include states which have Z-errors at both boundary 
sites, and in this case XE'^^^ — 2A. 

The second order correction is calculated in an iden- 
tical manner to the case with periodic boundary condi- 
tions. We have 



0(2) = 



PlVPiVPl 2NsPl + 2J2 



(A 



(0) 



E- 



-25 



(72) 



The above operator is already diagonal in our chosen ba- 
sis (the states of the form |C{a, (3, . . .})). Of the states in 
T, the cluster state |C) is the unique -1-1 eigenstate of all 
the stabilisers in the above sum, and so will be the lowest 
eigenvalue of 0'-^-'. Thus the second-order correction for 
the energy associated with this state is 



^ ^\C) - ^ 



-2Nsj . (73) 



So in the case of fixed boundary conditions, the cluster 
state is still the ground state produced (to zeroth order), 
with energy to 2nd order of 



E\c) = -gNs ~2X- 2Ns 



(74) 



A state |C{a, /3, . . .}) with ng Z-errors at boundary sites 
and nj Z-errors at interior sites will have energy: 

E{nB,n,) = E\c) + ns^B + ^^/A/ (75) 

where A/ := 2X^ / g and := 2A. Provided X/g < 1, 
the energy gap will be A/ — 2X^/g, the same energy gap 
which was obtained with the periodic boundary condi- 
tions. 

The first order corrections to this ground state will be 
given by 



A 



9 



(76) 



tJ.^BuB2 



and the bound on d, the average fidelity per site, remains 
the same as for periodic boundary conditions. 



B. Square lattice with fixed boundaries 

Consider a square lattice of dimension I x I. We define 
the number of interior (non-boundary) sites to be Ns, 



and so I = \/Ns + 2 and the total number of sites is 
Ns + 'iy/Ns + 4. Interior sites have coordination number 
4, edge sites coordination number 3, and corner sites co- 
ordination number 2, determining the number of qubits 
at each site. Denote the set of corner boundary sites 
by £1, the set of edge boundary sites by £2, and the 
set of interior sites by £3. Each of the three type of 
sites will have a different unperturbed spectrum at each, 
corresponding to Figs. [2l [3] and [4] respectively. The ze- 
roth order energies of the lattice are now much more 
complicated due to the presence of these three differ- 
ent types of sites. The ground state of gHo for the 
entire lattice, spanned by all possible logical states, is 
2'^^^+''^'^^+^)-dimensional ground-state space with en- 
ergy E^"^ = -4:gNs - I2gy/N^ - 8g. The next four ex- 
cited states separated by a energy gaps of 2g. 

At first-order in the perturbation, 6''-'^' — PlVPl = 
0, and thus there is still no first-order correction to the 
energies. At second, third, and fourth order, we have 



0(2) 



(jVs + 3V7Vs + 4)FL+E^g£,^A. 



0(3) = _ 



3!E 



(-4ff)2 



0(4) 



Wg-^ 



■ An 



49» 



1 
2? 



S^Si, 



15/16A^g+45/16v/7Vi^+10-(Ng/2+3^/iV^/2+4)nl p 



(77) 



(In this expression, m is the number of Z-errors at sites 
in £1 relative to the cluster state. It appears is this 
expression because O'^'^^ depends on the second-order en- 
ergies.) This operator is diagonal in the familiar basis 
{|C), |C{a}), \C{a, /?}), . . .}. The corresponding correc- 
tions to the lowest energy ground-state energy are 



x^e'^J^ 



1C) 



\c) - 2g2 
4 ^(4) _ -6iVs 



128 A" 



16 



5" 



(78) 



Thus, the non-degenerate ground state of the system is 
the cluster state |C), to zeroth order, with an energy to 
fourth-order given by 



So = - 4<?iVs - Ug^Ns - 8g - 

3v^,. 6Ns + 3VN^ 




^X\ 



(79) 



A state \C{a, /?, . . .}) obtained from the cluster state by 
ni Z-errors at sites in £1, 77.2 Z-errors at sites in £2 and 
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713 ^-errors at sites in £3 will have energy 



(ni,n2,ri3) 



E\c) + A(ni) + 712A2 + ngAg 



(80) 



where 

A3 
A(ni) := ( 



5A^ 

9 g- 

:= , ki 



A2 
6A4 



3A3 
452 

+ — 



-ni 4 

' g' 



(81) 



8A^ 

5^ 



For the range X/g < 1, the energy gap will be A3 — 
{5/8)X'^/g^, the same energy gap obtained using periodic 
boundary conditions. 

In summary, the perturbative procedure is still suc- 
cessful on the line and square lattice with fixed bound- 
ary conditions, producing an approximate cluster state 
as the non-degenerate ground state. Furthermore, the 
energy gap to the first-excited space is unchanged from 
the case with periodic boundary conditions. Extending 
these results to other lattices is straightforward. 



VI. THE CLUSTER STATE ON A GENERAL 
GRAPH 

We have shown that the perturbative procedure pre- 
sented here is successful in producing a non-degenerate 
ground state that approximates the cluster state on all 
of the lattice types examined so far. In fact, the cluster 
state on any lattice type with any boundary conditions, 
or more generally on any graph, can be approximated 
using this method. We now outline a proof of this result. 

For any graph, we place at each site a number of 
physical qubits equal to the coordination number (the 
number of bonds connecting that site to others) and 
take the Hamiltonian as in Eq. ([6|). We note that 
the form of the site Hamiltonian Hq needs only yield 
a two-dimensional degenerate ground state spanned by 
|00...0) and |ll...l) of all qubits at each site; aside 
from this requirement, its precise form is quite flexible. 

Wc first show that the operators produced at each or- 
der in perturbation theory will always possess the cluster 
state as an eigenstate, and more generally are diagonal- 
ized by the set of cluster states with Z errors. Note that 
the operators that arise at each order of the perturbation 
theory are linear combinations of terms of the form 



PLV[Ui=A^"''PLV)]PL, 



(82) 



for some integers I and ak, where Q :— (iJg"'' — gHg) ^. 
Now, y is a sum of operators cr^^ (g) crf^ j-^ where /i ^ i^, 

and therefore PLV{lli=ii^°"' PlV))Pl wiU be a sum of 
operators which map states in the logical space through 
the illogical spaces (by applications of cr^ over 
various bonds) and then return it to the logical space. 



From the Pauli operator commutation relations, we note 
that every term (g) that is applied to the logical 
space either commutes or anticommutes with the terms 
in Hq and so it will always yield an eigenstate of Hq. 
Thus, successive applications maps the logical space to 
eigenspaces of definite unperturbed energy, and the term 
fl = {Eq^^ —gHQ)~^ will just be a multiplicative constant. 

Furthermore, the fact that each application of cr^ (g) cr^ 
keeps the system in some eigenstate of Hq means that all 
the terms in the sum are of the form P^KPl, where K is 
some product of the . Now, suppose K does not 

commute with all the terms in Hq (i.e. it anticommutes 
with at least one of them), then the effect of this term 
will be that K maps logical states to illogical ones (the 
resultant state will have a —1 eigenvalue for at least one 
term in Hq, whereas the logical space is the -1-1 eigenstate 
of all the terms in Hq). Thus in this case PlKPl — 0. 
The only non-zero operators in the perturbation theory 
will be of the form P^KPl where K commutes with all 



But then. 



the terms in Hq (i.e. all the a 

if K commutes with each term in Hq it commutes with 
Pl, the projection onto the logical subspace. Thus, we 
have the eigenvalue relation 

PlKPl\C) = PlKPl [Pl\C2)\C2) ■ ■ ■ IC2)] 
= PlK\C2)\C2)---\C2) 
= Pl\C2)\C2)---\C2) ^\C) (83) 

where the last line follows because cr^ (g) cr^|C2) — ® 
f^|C2) = IC2), and K is a product of cr^ (g tr^ terms. 

Hence, we have shown that all the terms that arise at 
each order in the perturbation theory stabilise the cluster 
state. This certainly shows that, to zeroth order, the 
cluster state is one of the eigenstates selected out of the 
degeneracy by the perturbation. We now show that it is 
the non-degenerate ground state. 

Each term of the form P^KPl stabilises the cluster 
state, and the eigenvalues of P^KPl are restricted to 
±1. Therefore, all that must be checked is that the sign 
in front of P^KPl in the perturbation theory is negative 
to ensure that the cluster state is selected as the ground 
state. Suppose X is a product oim ® terms. Then 
the operator P^KPl will first appear at the mth order 
as a term in PlVIVIPlVY^-^Pl. The operators will 
always contribute a sign (—1)™^^ to the energy correc- 
tion. Furthermore, each carries with it a negative 
sign in the definition of V , contributing a further (—1)"* 
to the the energy correction. Therefore, each P^KPl 
will always appear in the energy correction with a neg- 
ative sign, thus selecting the cluster state as the ground 
state. Moreover, because it is clear that the cluster state 
stabilisers S"^ = (giy^^ Z,y, will always arise as one of 
the PlKPl terms in the perturbation theory, the ground 
state must be non-degenerate (as the state stabilised by 
these operators is unique). Hence we have shown what 
we set out to prove: that the cluster state on any lattice 
type can be approximated using this method. In fact, this 
perturbative approach can also be further generalised to 
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approximate other states with a PEPS description. 



VII. CONCLUSION 

The existence of gapped quantum many-body systems, 
with Hamihonians consisting of only two-body nearest- 
neighbour interactions, for which the ground state en- 
codes a chister state allowing universal MBQC is an ex- 
citing result for the potential realisation of a quantum 
computer. The obvious avenue for future investigation is 
whether existing natural or artificial materials exist with 
interactions similar to those described here. 

As we have shown, as the perturbation parameter \/ g 
becomes larger, the ground state begins to deviate from 
the cluster state due to perturbative corrections. In this 
work, we have analysed these corrections as a source of 
error. It may also be fruitful, however, to analyse the 
usefulness of the finite \/g ground state for MBQC in 
terms of the performance of a universal set of quantum 
gates, as in [24i|. Although we do not believe the model 
investigated here exhibits a phase transition at any \/ g 
(this is however an open question), it may nevertheless 
be possible that the usefulness of the ground state for 
MBQC undergoes some form of sharp transition [25| . 
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APPENDIX A: PERTURBATION THEORY 

We briefly outline the formalism of degenerate pertur- 
bation theory and the notation that we use, closely fol- 
lowing Ref. |22| . 

Suppose that the Hamiltonian has the form Hq + V, 
and that the eigenvalue problem has been solved exactly 
for Hq. The corrections brought about by the introduc- 
tion of the perturbation V can then be approximated by 
a power series expansion in V. For perturbation theory 
to converge, the magnitude of the largest eigenvalue of V 
must be smaller than that of Hq. 

Suppose the unperturbed spectrum has a degenerate 
subspace Ti.^ with energy E^^^ and we are interested in 
finding out how this energy degeneracy is broken. After 
the perturbation has been apphed, denote the perturbed 
eigenstates of this subspace by \^pi) and the perturbed 
energies hy Ei, for i = 1, . . . ,dim7ii, i.e.. 



Denote the projection onto the degenerate subspace Hl 
as Pl and define Pl = I — Pl- Then we can decompose 
liJi) as = |V'i)L + where ■= Pl^Ji) and 

Applying this decomposition to Eq. (|A1[) . we have: 

{Ho + V)\^P,)l + [Ho + V)\^,)i^ E,\,P,)l + E,\^,)i. 

(A2) 

Note that PlHo = HqPl = E^^ Pl and therefore that 
PlHq = HqPl. Multiplying Eq. jAl by Pl and Pl 
respectively, we obtain 

{E, - Ef^ - PLVPLm)L - iPLVPLM)L - (A3) 

{E, - Ho-PLVPLm)L- (PlVPlM)l = 0. 

(A4) 

Eq. (|A4p has the formal solution 

\^^)L - iE^ - Ho-PLVPL)-HPLVPLm)L , (A5) 

which can be substituted back into Eq. (|A3p to obtain 



where 



0|V'z>L = (i?.-i?f )I^.>L, (A6) 



e := PlVPl 

+ (PlVPlKE, -Ha- PlVPl)-\PlVPl) . (A7) 

This equation allows us to determine perturbed energy at 
any order of the perturbation theory. So far no approxi- 
mations have been made. To implement the perturbation 
theory, it is necessary to expand as a power series in 
V. We use 

iE,-Ho-PLVPLr^ 

= {I- {E, - Ho)-^PLVPLy\E, - Ho)-^ 

OO 

= ij2 - HQr^PLVPLy"){E. - Ho)-^ ■ 



m=0 



(A8) 



The energies Ei in these expressions must also be ex- 
panded as power series, Ei = Ej^ 
we have 



(0) 



EZiEl'l Then 



{Ei—Ho 



{eT 



{k) 



Hn 



(k) 



{Ho + V)\ij,) =E,\iJ;, 



(Al) 



= /+(4"'-Ho)-^^K 

k=l 

= An 

where we have defined the operators 
n := [E^p - Ho)-' , 

OO 

A:= (/ + f7^i?f)) 



(4"' - Ho)-' 
(A9) 

(AlO) 
(All) 



fc=i 
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The operator 9 can then be expressed as 

oo 

9 = PlVPl + PlV [^^PLVr^^PLVPL , (A12) 

m— 

where we have used the fact that P^ commutes with Q, 
and A. Additionally, A must be further expanded out as 
a power series 

oo oo oo 

A = (/ + f7 ^ E^r' = Y,{-n ^ E['^y . (A13) 

k=l j=0 k=l 

We can now identify the terms in Eq. (|A12[) of each 
order. Specifically, denote the terms in 9 of fcth order by 
6'('^' , so that 9 = X^fc^i • Then, when approximated 
to nth order in V, Eq. (jA6p becomes 

n n 

^0W|^,)^ = ^£;W|^,)^, (A14) 
fc=i fc=i 

which is an eigenvalue equation over the subspace L. (We 
note that \tpi)L 7^ because \4>i) G H-l at zeroth order 
and the perturbation is assumed small). The energy cor- 
rections to nth order -^1*^^ eigenvalues of the 
operator X)fc=i ^^'^^ ■ Furthermore, the eigenstates that 
are selected by the perturbation to break the degeneracy 
(to zeroth order) are just the eigenvectors of J2k=i ^^^'^ 
corresponding to each eigenvalue. Note that 6*^*^^ depends 
on lower order energy corrections and so each lower or- 
der correction must be calculated before proceeding to 
the higher order corrections. At each stage we must in- 
sist that the nth order energies differ from the (n— l)th 
order energies only by an amount nth order in V , which 
removes any non-physical solutions. 

To determine the explicit form of the 9'^''\ we sim- 
ply substitute Eq. (|A13p into Eq. (|A12|) and identify the 
terms of the required order. Clearly 9'^'^^ = PlVPl, and 
so to first order in V, Eq. (|A14[) becomes 

0^''>\^^)L = PlVPl\^,)l = . (A15) 

Thus, the first order energy corrections to states in Ti.^ 
are just the eigenvalues of the matrix PlVPl- In the 
non-degenerate case we see that Eq. (|A15p reduces to the 
well-known expression £:f) = (^l^^l^)- If this eigenvalue 
spectrum is still degenerate, then the degeneracy is not 
completely broken at first order. It is then necessary to 
go to second order 

{E^'^ +eP)\^,)l^{9^'^ +9^^^)\^,)l 

= (PlVPl + PlVVlPlVPlMCjl ■ 

(A16) 

Again we can examine whether the degeneracy has 
been broken at this stage by examining the eigenval- 
ues and eigenvectors of the above operator PlVPl + 
PlVUPlVPl- If not, we must continue to proceed to 



higher orders. The formulae for higher orders become in- 
creasingly complex, but they simplify if we assume that 
e'^^^ — 0, which will be true for all the cases of interest 
that we investigate in this paper. For example we find 
(with ^ = 0) 

0(1) = PlVPl 

6»(2) = PlVOPlVPl 

6»(3) = PLV[nPLV]^PL 

6»(4) = PLV[nPLV]^PL - E^'^^PlVQ^ PlVPl , (A17) 
and so forth. 

There are two additional points that must be noted. 
First, from Eq. (jAl2|l . it can be concluded that 

J2k=i ^i''^ the eigenvalues of the operator X^fcLi 

(k) 

for each order N. However, the energies El are not 
generally the eigenvalues of the operators 6*''^' ; this will 
only be true in general when all operators 0^'^) can be si- 
multaneously diagonalised. Fortunately, for the systems 
investigated in this paper, it can be shown that the ^^'^^ 

Ik) 

commute with each other and therefore the energies E^ 
are indeed the eigenvalues of the operators 6''-'^-' . Second, 
there are some general properties we can note about the 
form of the 0^'^) operators. From the form of Eq. ()A12|) 
and the series expansion in Eq. (jA13P , it is clear that all 
terms in the expansion of 9 are proportional to operators 

of the form PL^^(nLi(^"'^i^))^i some l,ak G N. 
This result is used in Sec. IVII 

The above analysis allows the zeroth order energy 
eigenstates to be determined. We now direct our atten- 
tion to the first order corrections to these states. Suppose 
we have determined that |V'i) = IV'i) l = \i) to zeroth or- 
der for some \i) G 1-Ll using the above method. From 
Eq. ((X5|) and using Eqs. and (IX9|) we have, to first 

order in V ^ 

mL = ^{PLVPLm)L= E J'^'L b')^ 

(A18) 

where \j) G Ti^ are the eigenstates of Hq with energy 
E^l^\ Determining \ipi)L to first order is somewhat more 
complicated: even though \4>i)L = \i) to zeroth order, it 
is possible for first order corrections to come from other 
states \l) G T-Ll- However, if the energy degeneracy be- 
tween \i) and |/) is only broken at order to/, then one 
must go to the order (m; -I- 1) equations just to deter- 
mine the first order eigenstate corrections io \i}}i)L- In 

this case, Eri7' = Y.k=i^ E^'^^ but e'^""^ ^ eI""\ 
Then, at order (toj -|- 1), Eq. (|A14p reads 

mi m/ + l 

(^eW +0("'+l))|^,)i = J2 4'^\^^)L■ (A19) 
k=l k=l 
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Taking the inner product with \l) and rearranging gives 

{i\e("^'+'^\A)L-El"''+'\i\A)L 



'"1 + 1)1 

(A20) 



to first order in V. Therefore, to first order in we have 
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